On the Physics of Size Selectivity 
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We demonstrate that two mechanisms used by biological ion channels to select particles by size 
are driven by entropy. With uncharged particles in an infinite cylinder, we show that a channel that 
attracts particles is small-particle selective and that a channel that repels water from the wall is 
large-particle selective. Comparing against extensive density-functional theory calculations of our 
model, we find that the main physics can be understood with surprisingly simple bulk models that 
neglect the confining geometry of the channel completely. 
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Ion channels are membrane-spanning proteins that 
passively transport ions down their electro-chemical gra- 
dients. Channels can open and close their pore upon 
stimulation, a process called gating. In addition, some 
channel types can preferentially select the ion species 
they conduct. These properties make channels respon- 
sible for a large number of physiological phenomena in- 
cluding propagating action potentials along neurons and 
initiating muscle contraction. Experimentally, channels 
can be studied either in vivo or in well controlled, but 
non-physiological conditions one channel at a time [ij . 

Experimentally, eleven "selectivity sequences" have 
been found that rank the relative preference of a channel 
for conducting alkali metal ions through cation channels 
In this paper we focus on the two extreme selectiv- 
ity sequences where the channel prefers to conduct small 
ions or large ions. Examples of small-alkali metal cation 
selective channels include the L-type calcium channel, 
the ryanodine receptor (RyR), and the neuronal sodium 
channel. Examples of large-ion selective channels include 
voltage-gated potassium channels, gramicidin A (gA), 
and the nicotinic acetyl choline receptor (nAchR). 

All channel pores are several nanometers in length and 
have radii ranging from approximately 2 A for gA and 
potassium channels to > 5 A for porin channels and 
nAchR. Although its dimensions are small and similar 
to the size of ions, equilibrium bulk models for selectiv- 
ity [^IMS that neglect the confining due to the pore can 
help to understand experiments and reproduces results 
from Monte Carlo simulations qualitatively [3- These 
and other studies of selectivity 0j 0, S found that 
highly-charged channels such as the sodium, L-type cal- 
cium, and RyR channels are small-ion selective because 
smaller ions balance the charge on the protein in a smaller 
volume, a mechanism named charge/space competition. 
For large-ion selectivity there are several theories, some 
based the dehydration energy of ions 0, 0| and others 
based on the repulsion of water from hydrophobic regions 
of the channel, without explicit ion hydration Q. 

In the present Letter we wish to extend the under- 



standing of the physical mechanism that underlies size 
selectivity by highlighting entropy as driving force in 
both selectivity mechanisms. To do so we choose a model 
with the smallest number of parameters that allows one 
to study the phenomenon of size selectivity. This model 
consists of a mixture of uncharged hard spheres with dif- 
ferent radii 11]. This model obviously neglects electro- 
statics as well as effects caused by hydration shells. In 
our model the hard-sphere "cations" are attracted into 
and the hard-sphere "anions" are repelled from the selec- 
tivity filter through an effective potential, which param- 
eterizes the long-ranged contribution such as the Don- 
nan potential. We denote its amplitude by Uattr > 0. 
We have verified that, due to the mentioned repulsion, 
the presence of anions in the system has a very small 
effect on the results. For the sake of simplicity we there- 
fore include only the solvent and hard-sphere "cations" 
in our model. For studying large-ion selective channels 
such as the gA channel we introduce an effective repul- 
sion Vrep > 0, which models the hydrophobic repulsion 
of water from the channel. 

Here we focus on the importance of the entropy of both 
the ions and the solvent, which is often neglected models 
of biological systems. Our approach is in line with the 
findings of theoretical studies of biological problems in 
different areas such as protein folding . If the solvent 
is modeled as a fluid of particles with a size comparable 
to that of the other species, it provides a crowded envi- 
ronment for the ions, which in turn have to compete for 
free space even at low ion concentrations. Uattr and Vrep, 
the amplitudes of the ion attraction and the water repul- 
sion, respectively, are simply parameters in our model. 
Since we wish to understand the phenomenon of size se- 
lectivity among equally charged ions, a detailed account 
of electrostatic and hydrophobic interactions will modify 
the selectivity of the channel quantitatively, but not the 
role of entropic forces described here. 

In order to test if such a simple system selects particles 
by size, we start by considering two compartments. The 
bath, denoted compartment /, contains ions at given con- 
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centrations of the order of lOOmM dissolved in a crowded 
solvent of 55.5 M, the density of pure water under nor- 
mal conditions. A second compartment, denoted com- 
partment //, models either the selectivity filter of the 
channel that attracts ions by the action of Uattr or a hy- 
drophobic channel that repels water by the action of Vrep- 
Inhomogeneities caused by the confining geometry of the 
protein are ignored at this stage, but are included later 
in our density functional theory (DFT) calculations. 

By allowing equilibrium between the compartments, 
the concentration of all components in the filter adjust 
so that the grand potential of the system is minimized. 
This is described by the equality of the chemical potential 
/if = df{{p"})/dpi of species i between compartments 
a = /, //, where is the corresponding number density 
of this component. For the free energy density of the 
mixture / we employ the expression correspondin g to 
the Mansoori Carnahan Starling and Leland (MCSL) |l3| 
equation of state. For the ion components we obtain 

^^l{{p]})=^^V{{p]'})+UaUr, (1) 

and for the water component 

^^ko{{p]}) = o({pf }) - Vr.,. (2) 

Uattr and Vrep are the electrostatic and hydrophobic 
parts of the chemical potential. Equations and |5J) are 
coupled, highly non-linear equations that determine the 
densities in the filter compartment pl^ for given densities 
in the bath, , and parameters Uattr and Kep- This bulk 
model allows us to study the physics of size selectivity 
assuming that inhomogeneities in the pore are unimpor- 
tant. In order to quantify which component is preferen- 
tially absorbed in the pore, we define the absorbance of 
component i, = pj^ /pf, which compares the density 
of component i in the filter to that in the bath. The 
selectivity Sij of the filter is defined by Sij = ^i/^j- 

The assumption that the inhomogeneities are not im- 
portant can be tested by taking the confining geometry 
of the pore into account within the framework of DFT. 
Here we employ a recently improved version of Rosen- 
feld's fundamental measure theory which is based on the 
accurate MCSL equation of state In addition to the 
inputs to our bulk model, the DFT approach also requires 
a model for the protein that forms the pore of the selec- 
tivity filter. The simplest way to incorporate the effect 
of the protein in the DFT approach would be to define 
external potentials Vf^'^{r) that confine particles of all 
components inside the channel. Unfortunately, the ac- 
tual form of Vj*^^* (r) is unknown because of the unknown 
structure and mechanical properties of the protein. In 
our study we choose to model the protein as a hard- 
sphere fluid which is restricted by a hard-wall potential 
to a region outside the pore. If the "protein fluid" is suf- 
ficiently dense, spheres of other components will find it 
difficult to enter the protein. In this way the protein fluid 




FIG. 1: The absorbance of Na and K in an attractive filter 
as function of Uattr- We compare bulk results (lines) for a 
binary mixture to those for a ternary mixture. In the ternary 
mixture of water, Na, and K the competition for space leads 
to a selectivity of the smaller component, SNa,K S> 1, as is 
shown by the solid line in the inset. The symbols, which are 
explained in the text, denote DFT results for the ternary case. 

acts as a rough wall that can be penetrated by particles 
of other components, which is more appropriate than a 
smooth, hard wall. We model the pore as an inflnitely 
long cylinder, which simplifles the calculations because 
the density proflles Pi{r) — pi{r) of all components, la- 
beled by i, depend only on the radial distance r. 

By minimizing the density functional of a N- 
component mixture J7[{p,i(r)}] we obtain the inhomoge- 
neous equilibrium density profiles Pi(r), as well as the 
grand potential fl of the system. From these we derive all 
quantities of interest. As mentioned above, we consider 
an external potential acting only on the protein compo- 
nent, which in turn acts on the other components in the 
pore so that they develop inhomogeneous structures. If 
we take the limit of the protein fluid density going to 
zero, no external potential acts on the system and we re- 
cover our bulk approach. In addition to the parameters 
entering our bulk approach, two new parameters are re- 
quired to describe the system fully, namely the protein 
fluid packing fraction rjp and the pore radius Rpore which 
defines the region r < Rpore from which the protein fluid 
is expelled. Here we chose rjp — 0.45 and the diameter of 
the particles that constitute the protein to be 2.45A. 

In our DFT calculations we obtain the absorbance of 
species i inside the pore from its density profile pi (r) by 

= ~rH5 — / drrpi{r). (3) 

Pi ^pore Jo 

If the inhomogeneities in pi{r) are small this definition 
yields numbers for the absorbance and the selectivity 
very close to those predicted by the bulk approach. 
We start by reporting results of our bulk approach 
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FIG. 2: Typical density distributions of Na (a) and K (b) in- 
side an attractive, cylindrical filter as obtained from DFT for 
various values of J7aiir . r denotes the radial coordinate. While 
the density of Na (a) increases monotonically with Uattr, the 
density of K (b) shows a non-monotonic behavior. For the 
radius of water we use Rh^o = <yH20 /2 — 1.4A 

which highhght the role of competition between ion 
species for free volume in the filter. We consider a bath 
consisting of a binary mixture of lOOmM Na or lOOmM 
K in water and attract ions into the filter by the poten- 
tial Uattr: the range from to 10 ksT. As a result 
of the attraction, the concentration of Na or K in the 
filter increases. We quantify the effect of the attractive 
potential by the absorbance ^i, i = Na,K, and show 
its dependence on Uattr in Fig. ^ The dotted line in 
Fig. n shows the result for the binary mixture of Na and 
water, obtained by the bulk approach, and the dashed- 
dotted line shows the corresponding quantity for the bi- 
nary mixture of K and water. Comparing these results, 
we observe a similar increase in the filter compartment 
for both species, especially at small values of Uattr, be- 
cause in that regime the particles still find free space 
between the solvent particles. At higher values of Uattr, 
however, Na must squeeze out less water from the fil- 
ter than K to follow the attraction, because a^a < f/f • 
The result is that S^Na > and hence the selectiv- 
ity SNa,K = ^Na/£,K > 1, SiS IS showu in thc inset of 
Fig. n (dotted line). An attractive channel favors smaller 
species. In this case there is no direct competition be- 
tween species of different size and the selectivity SNa,K 
remains small, leveling off at 2.2 for large values of Uattr- 
The situation becomes much more interesting and 
richer if we consider a ternary mixture of water, lOOmM 
Na, and lOOmM K in the bath. Both ions species feel the 
attraction and are in direct competition for the free vol- 
ume in the filter. The result of this competition, within 
the bulk approach, can be seen in Fig. ^ (full and dashed 
lines). For weak attractions the situation is similar to 
the previous system, i.e. the ions fill the free volume 
between solvent particles and the actual size of the ions 
is unimportant. At stronger attractions, however, when 



the ions no longer find free volume, they have to squeeze 
water out of the filter to follow the attraction. When the 
water concentration in the filter is significantly reduced, 
direct competition between the ion species sets in and 
it becomes energetically and entropically more favorable 
for the system to allow the smaller ion species to fill the 
volume by squeezing out the larger ions. The Na con- 
centration in the filter continues to increase while the K 
concentration starts to decrease. If the attraction is suf- 
ficiently strong, this effect can be dramatic and results 
in a large selectivity SNa,K ^ 1, as we show in the in- 
set of Fig. n (full line) . It is important to realize that 
this physical mechanism always prefers the smallest ion 
species and prefers Li over Na if we were to start with a 
bath consisting of water, Li, and Na. 

We test the predictions of our bulk approach by ex- 
tensive DFT calculations of our model by changing the 
pore radius Rpore from 1 to SA. Some DFT results for 
-Rpore=2.1 A (square), 3.5 A (circle) and 5 A (diamond) 
are shown in Fig. ^ and its inset. Full and open sym- 
bols denote results for K and Na, respectively. We find 
that if the pore is sufficiently wide, the DFT results for 
the selectivity are in good agreement with our bulk pre- 
dictions, although the absorbances obtained from DFT 
deviate slightly from the corresponding bulk values. This 
agreement indicates that in this case the inhomogeneities 
caused by the confinement are moderate. Most chan- 
nels whose radii are known have radii >3 A P, 0, 0|. 
If, however, the pore radius becomes smaller (< 2.5 A) 
and nearly equal to the particle radii, the confinement 
becomes increasingly important and results in deviations 
between DFT and bulk results. In this regime, we observe 
nonlinear absorbance behavior similar to that described 
by Goulding et al. In narrow attractive channels 

the small ion selectivity can be increased dramatically. 
Channels known to have radii this small include gA [l| 
and the KcsA potassium channel 1^. 

A set of density profiles of the Na and K species for 
various values of Uattr is shown in Fig. El Thc pore radius 
is Rpore = 3.5A and the ion concentrations in the bath 
are pNa = Pk = lOOmM. These profiles demonstrate the 
monotonic increase of Na and the non-monotonic behav- 
ior of K in the pore as Uattr is increased. We find that 
the structure in the density distribution of the ions follow 
closely the densities predicted by the bulk approach. 

Now we turn to a different type of channel that is 
characterized by hydrophobic protein walls, which we 
take into account by the effective water repulsion Vrep- 
Weakly-charged channels such as gA and nAchR seem 
to have these properties 0, and references therein]. The 
little charge in these channels is sufficient to repel anions 
from the channel but not to distinguish between cations 
of different size. In order to keep our model as simple 
as possible, we neglect the anions and set Uattr = 0, 
however, we have verified that an extended model that 
includes a weak cation attraction and hard-sphere "an- 
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FIG. 3: The absorbance of Cs and Na in a hydropho- 
bic fiher as function of Vrep- With increasing value of Vrep 
both ion densities increase, however, the density of the larger 
component increases faster. The selectivity of Cs over Na is 
shown in the inset. Lines denote bulk results, while symbols, 
which are explained in the text, those obtained by DPT. 



ions" which are repelled from the channel predict equiva- 
lent results. Following Ref. we consider Vrep < 'iksT. 
In the present study we compare lOOmM of both Na and 
Cs, which have a more pronounced size difference than 
Na and K considered earlier. 

In Fig. Owe show the absorbance of Na (full line) and 
Cs (dashed line) in a hydrophobic pore as function of 
Vrep for a ternary mixture with pi^a = PCs = lOOmM, 
predicted by the bulk theory (lines) and the correspond- 
ing selectivity Scs,Na in the inset. While the densities of 
both species increase monotonically as the water repul- 
sion increases, the density of the larger component, Cs, 
increases faster than that of the smaller component, Na. 
Hence the channel is large ion selective. 

Again we test the prediction of the bulk approach by 
DFT calculations. In order to account for the hydropho- 
bic interaction between the protein and water we intro- 
duce an external potential acting on the water. V^'!^Q{r) 
repels water only from a 1.4 A neighborhood of the pro- 
tein. This is a clear difference between the DFT and the 
bulk approach, which by construction cannot take surface 
effects into account. This difference is reflected by lower 
Cs densities in the pore obtained from DFT (symbols) 
as compared to the bulk theory (lines) in Fig. O and the 
smaller selectivity in the inset of Fig.|21 Here the squares, 
circles and diamonds denote DFT results for Rpore =2.8, 
3.5 and 5 A, respectively. This difference gets smaller as 
the radius of the pore is reduced. The overall agreement 
agreement is still qualitatively good. 



Our extensive DFT study confirms the validity of the 
bulk approach for both kinds of channels also for different 
choices for rjp and the ion concentrations. The confining 
geometry of the pore becomes important only if the ra- 
dius of the pore is comparable to the radii of particles 
inside the pore. Our findings also show clearly that the 
entropy of the mixture of particles is enough to give size 
selectivity and hence has to be taken into account prop- 
erly. A very important, yet often neglected, contribution 
to this entropy stems from the fact that water is a dense 
fluid and leaves little space to the ions. 
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